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ABSTRACT 

Variable-delay Polarization Modulators (VPMs) are currently being implemented in experiments 
designed to measure the polarization of the cosmic microwave background on large angular scales 
because of their capability for providing rapid, front-end polarization modulation and control over 
systematic errors. Despite the advantages provided by the VPM, it is important to identify and 
mitigate any time-varying effects that leak into the synchronously modulated component of the signal. 

In this paper, the effect of emission from a 300 K VPM on the system performance is considered and 
addressed. Though instrument design can greatly reduce the influence of modulated VPM emission, 
some residual modulated signal is expected. VPM emission is treated in the presence of rotational 
misalignments and temperature variation. Simulations of time-ordered data are used to evaluate the 
effect of these residual errors on the power spectrum. The analysis and modeling in this paper guides 
experimentalists on the critical aspects of observations using VPMs as front-end modulators. By 
implementing the characterizations and controls as described, front-end VPM modulation can be very 
powerful for mitigating 1// noise in large angular scale polarimetric surveys. None of the systematic 
errors studied fundamentally limit the detection and characterization of B-modes on large scales for 
a tensor-to-scalar ratio of r = 0.01. Indeed, r < 0.01 is achievable with commensurately improved 
characterizations and controls. 


1. INTRODUCTION 

Cosmic microwave background (CMB) measurements 
have been the cornerstone of cosmological advances. This 
is particularly true for measurements that test and con¬ 
strain inflation. The COBE space mission discovered the 
anisotropy and measured the amplitude of scalar fluctu¬ 
ations (Smoot et al. 1992; Bennett et al. 1994). The 
WMAP space mission confirmed the COBE scalar am¬ 
plitude and reported the predicted small deviation from 
scale invariance with the scalar spectral index, < 1, 
and flatness to within a fraction of a percent (Hinshaw 
et al. 2013; Bennett et al. 2013). The Planck space mis¬ 
sion confirmed the WMAP results and placed tight new 
limits on the gaussianity of the fluctuations (Planck Col¬ 
laboration et al. 2015). Looking forward, the next chal¬ 
lenging step is to test and constrain inflation by charac¬ 
terizing the predicted stochastic background of gravita¬ 
tional waves via their imprint onto the 5-mode polariza¬ 
tion of the CMB. Measurement of this signal would es¬ 
tablish the energy scale of inflation. It remains to be seen 
if characterization can be done from either the ground or 
a balloon, especially at large angular scales, or if it will 
require another space mission. The current best 95% up¬ 
per limit on r, the tensor-to-scalar ratio, is r < 0.07 at 
k = 0.05 Mpc“^ (BICEP2 and Keck Array Collabora¬ 
tions et al. 2016). 

On angular scales less than ^ 2 ° {£ > 100), the cos¬ 
mological 5-mode angular power spectrum is dominated 
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by the gravitational lensing of the curl-free 5-mode po¬ 
larization spectrum that is generated by density per¬ 
turbations in the early universe (Zaldarriaga & Seljak 
1998; Hu & Okamoto 2002). At larger angular scales, 
5-modes from inflationary gravitational waves are pre¬ 
dicted to source two distinct features in the angular 
power spectrum: an enhancement at ^ ^ 100 correspond¬ 
ing to the horizon size at recombination and another en¬ 
hancement at low < 10 corresponding to the size of 
the horizon at reionization (Zaldarriaga & Seljak 1997; 
Kamionkowski et al. 1997). A compelling detection of a 
gravitational wave background requires the characteriza¬ 
tion of 5-modes at angular scales £ < 200. Such a de¬ 
tection would not only test the amplitude of the posited 
signal, but also the predicted shape of the 5-mode an¬ 
gular power spectrum. 

Measuring the polarization is challenging due to the 
minute signal level and the fact that this signature must 
be characterized in the presence of relatively large back¬ 
grounds. These backgrounds include the unpolarized 
radiation from the CMB, Galactic emission, emission 
from the instrument, and in the case of ground-based 
and balloon-borne instruments, emission from the atmo¬ 
sphere and the ground. A key aspect of CMB polarime- 
try is the use of temporal, spatial, or phase modulation 
to enable separation of the desired polarized signal from 
residual environmental influences on the instrumental re¬ 
sponse. 

Modulation is the process of encoding signals of inter¬ 
est in a manner they can be subsequently synchronously 
demodulated and detected with high fidelity. Instru¬ 
ment influences outside the signal band and reference 
phase can be highly suppressed by the demodulation pro¬ 
cess. The process of modulation is inherently a “lock-in” 
measurement, and its use in microwave radiometry dates 


2 


back to Dicke (1946). In this signal processing technique, 
a modest degradation in sensitivity is traded for a dra¬ 
matically increased signal stability and/or control over 
environmental effects. 

A well-designed instrument and observing strategy 
commonly employs multiple levels of modulation at dif¬ 
fering frequency to aid in the separation of the desired 
signal from spurious artifacts. These levels are distin¬ 
guished by their characteristic time scales. Generally, the 
primary modulation, the one implemented on the short¬ 
est time scales, is critical because it sets the instrument’s 
ability to mitigate variations in the instrument and en¬ 
vironment. Such variations tend to be more severe at 
low frequencies and are typically characterized by a 1 //^ 
functional form, where a is the spectral index of the vari¬ 
ations. The primary modulation frequency is chosen to 
be above the knee frequency of the 1//^ spectrum. By 
only measuring the signal at the modulation frequency 
and phase (by demodulating and detecting it), the in¬ 
strument limits noise contributions outside of this band 
and specifically suppresses the effect of the variations in 
the system and environment that occur at lower frequen¬ 
cies than that of the modulation. Secondary modulation 
strategies are generally employed to remove systematic 
effects from the targeted signal. For example, instru¬ 
ment rotation can be employed to separate an astronom¬ 
ical polarization that is fixed in sky coordinates from an 
instrumental one that is constant in instrument coordi¬ 
nates (Keating et al. 2003; Hinderks et al. 2009; Barkats 
et al. 2014; BICEP2 Collaboration et al. 2014a; BICEP2 
and Keck Array Collaborations et al. 2015). 

The method of primary modulation most commonly 
utilized in CMB polarization experiments is an evolution 
of that developed in the context of continuum radiometry 
and photometry (Emerson et al. 1979; Emerson & Graeve 
1988; Mangum et al. 2007). These approaches have been 
successully applied in mapping the CMB anisotropy - 
the telescope is scanned in angle over a patch on the 
sky such that spatial frequencies are modulated and sub¬ 
sequently encoded at temporal frequencies in the data 
(Wright et al. 1996; Tegmark 1997; Pryke et al. 2009; 
Story et al. 2013; Diinner et al. 2013; Bennett et al. 2013; 
BICEP2 Collaboration et al. 2014b; POLARBEAR Col¬ 
laboration et al. 2014). Limits on telescope capabilities 
and variations in the unpolarized emission from the sky 
over large angular scales have typically limited suborbital 
scanning experiments to multipoles corresponding to the 
recombination feature and higher. 

An alternative to spatial scan modulation is to di¬ 
rectly encode the incident polarization signal. System 
level implementations of polarization modulation employ 
such devices such as half-wave plates (HWP) (Reichborn- 
Kjennerud et al. 2010; Eraisse et al. 2013; Kusaka et al. 
2014), Earaday rotators (Moyerman et al. 2013), and 
phase switched receiver architectures (Jarosik et al. 2003; 
Bersanelli et al. 2010; QUIET Collaboration et al. 2013) 
that vary the polarization state of the incoming radia¬ 
tion for synchronous demodulation and detection of the 
signal. In the ideal case, unpolarized and polarized com¬ 
ponents of the incoming radiation are treated indepen¬ 
dently by the system architecture and analysis pipeline. 
This approach enables configurations that approximate 
a null result when observing an unpolarized source. This 
is of particular value for mitigation of the effectively un¬ 


polarized atmospheric emission on large angular scales. 
However, the utility of a given modulation technique ul¬ 
timately depends on how well one can limit non-ideal 
performance during observations and calibration and, if 
necessary, model and remove such behavior via an unbi¬ 
ased procedure to achieve the desired measurement pre¬ 
cision. 

VPMs operate by introducing a variable phase delay 
between two orthogonal linear polarizations (see Eigure 2 
and Chuss et al. (2006)). Such modulators are imple¬ 
mented by placing a polarizing grid in front of and paral¬ 
lel to a flat mirror. The phase, and therefore the polariza¬ 
tion sensitivity of the instrument, is controlled by varying 
the separation between the grid and the mirror. VPMs 
are being applied to the measurements of the CMB polar¬ 
ization at large scales. This has many advantages. Eirst, 
VPMs can be made large enough to be usable as the first 
optical element for meter-sized telescopes. This mini¬ 
mizes systematic errors due to temperature to polariza¬ 
tion leakage in the instrument. Second, for CMB wave¬ 
lengths and meter-sized apertures, the modulation can 
be physically implemented at a frequency of several Hz, 
which moves the signal above the typical knee frequency 
of the sky noise for ground-based observations. Addition¬ 
ally, the VPM allows independent measurements of po¬ 
larization at each beam position without relying on scan¬ 
ning to reject the excess low-frequency noise. East polar¬ 
ization modulation can enable map-making without first 
differencing time-ordered data (TOD) from co-pointed 
orthogonally polarized detectors to remove the unpolar¬ 
ized signal. VPMs, due to the fact that they modulate 
polarization with linear-circular conversion, reduce sig¬ 
nal contamination due to cross-polarization {Q ^ U) 
compared to wave plates (Chuss et al. 2010). Since the 
circular polarization in the CMB is expected to be low 
(Kosowsky 1996), residual mixing between the linear and 
circular terms in the VPM transfer function will not 
cause significant systematic contamination of the linearly 
polarized CMB maps. 

The advantages of a VPM for recovery of polarization 
at large angular scales in the presence of atmospheric 
variability are captured in the maps of Eigure 1, which 
are generated following the simulation techniques dis¬ 
cussed in Section 4. The only difference between the 
simulations resulting in the two plots is the presence of 
a VPM, which modulates the sky signal before the ad¬ 
dition of the polarized 1// noise. Without a VPM, the 
1// noise contaminates the reconstructed map. With a 
VPM, the modulation allows the sky signal to be recov¬ 
ered on large angular scales. 

Ideally, a polarimeter would return a null demodulated 
signal or “offset” when viewing an unpolarized source. 
In reality, a measurement system has a finite instrumen¬ 
tal and environmental contributions to response which 
in practice one strives to minimize their spatial-temporal 
variation. These “systematic errors” must be character¬ 
ized and removed in order to meaningfully interpret the 
data. This is especially important in CMB polarization 
measurements for which the signal of interest is small, 
extended, and measured in the presence of large back¬ 
grounds. Detailed studies of the instrumental response 
have been previously carried out for a variety of sys¬ 
tem architectures (Shimon et al. 2008; MacTavish et al. 
2008; Miller et al. 2009; Karakci et al. 2013; Kusaka et al. 
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Fig. 1. — Recovered map of Stokes Q for simulations with polarized 1// noise. The simulations and reconstruction follow the description 
in Section 4. The amplitude of the noise has been arbitrarily chosen to illustrate the effectiveness of polarization modulation for large scale 
maps. A detailed study that employs a more realistic noise level is treated in Section 4. The leftmost panel shows the input Stokes Q map. 


2014). This paper undertakes a similar analysis for a 
VPM-based polarimeter architecture. In this work we 
focus on systematic influences arising from the radiative 
properties of the modulator, and develop techniques for 
separating these instrumental artifacts from the desired 
astrophysical signal. 

VPMs will be deployed on two upcoming CMB ex¬ 
periments. The first experiment is the Cosmology 
Large Angular Scale Surveyor (CLASS). CLASS is cur¬ 
rently being deployed to the Atacama desert in Chile 
(Essinger-Hileman et al. 2014). CLASS will consist of 
one 38 GHz telescope, two 93 GHz telescopes, and one 
dichroic 148/217 GHz telescope. Each telescope will have 
a VPM at ^ 280 K as its first optical element and they 
will observe ^ 70% of the sky from their —23° (south 
latitude) location in Chile. VPMs at 1.5 K will be de¬ 
ployed on the Primordial Inflation Polarization Explorer 
(PIPER). PIPER is a balloon experiment (Lazear et al. 
2014) that will operate at 200, 270, 350, and 600 GHz. 

This work uses simulations to propagate a realistic 
model of the instrument’s response in the time domain 
through polarized power spectra. The structure of this 
paper is as follows. Section 2 gives an overview of the ef¬ 
fect of the VPM on the input signal. Section 3 introduces 
systematic effects generated by the VPM. Section 4 de¬ 
scribes the simulations we perform. Section 5 shows the 
results of those simulations. Einally, we give the conclu¬ 
sions in Section 6. 

2. VPM MODULATION 

The VPM modulates polarization by introducing a 
controlled, variable phase delay between orthogonal lin¬ 
ear polarization states. As shown in Eigure 2, a common 
implementation of a VPM consists of a wire grid polar¬ 
izer placed in front of and parallel to a translating mirror. 
Light polarized parallel to the grid wires reflects off of 
the polarizer; light that is polarized perpendicular to the 
wires is transmitted through the grid and reflected off of 
the mirror. The two polarizations are recombined at the 
output with a relative phase delay that is proportional 
to the grid-mirror physical separation, z. 

One can visualize the effect of this introduced phase 
delay as follows. When the phase delay is zero, the out¬ 
put polarization is unchanged. As the phase delay is 
increased to 7r/2, the VPM becomes a reflective quarter- 
wave plate, converting linear polarization to circular po¬ 
larization. Increasing the phase delay further to tt, one 
recovers the behavior of a reflective HWP for which the 
direction of linear polarization is reflected about the line 
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Fig. 2.— A diagram of the VPM. Light that is polarized parallel 
to the grid wires is reflected. Light that is polarized perpendicular 
to the wires is transmitted through the grid and reflected off of 
the mirror. The two polarizations are recombined at the output 
with a relative phase delay that is a function of the grid-mirror 
physical separation, z. By modulating 2 ;, the polarization state is 
modulated. For CMB frequencies, 2 : will range from ps 0.5 — 5 mm. 


parallel to the wires. 

2.1. An Ideal VPM 

The key characteristic of an ideal VPM is that each 
polarization is treated by the VPM independently with 
100% transmission for the light polarized perpendicular 
to the wires and 100% reflection for the light polarized 
parallel to the wires. Thus, the ideal system is lossless 
with perfect polarization isolation. Defining a coordinate 
system such that the detectors are sensitive to Stokes Q, 
we assume an architecture in which the grid wires are 
at ±45° with respect to the orientation of the detectors’ 
polarization sensitivity. The polarization transformation 
of the VPM for monochromatic light can be written as 

s' = Ms (1) 

where M is the Mueller matrix and s and s' are the 
vectors of Stokes parameters before and after the trans¬ 
formation, respectively. Writing out the vectors and ma- 
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trices gives us: 
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Here, (j) = (j){z) is the phase delay and 2 : is the grid- 
mirror separation. The primed and unprimed Stokes pa¬ 
rameters represent the polarization states on the input 
and the output side of the VPM, respectively. For the 
case in which the wavelength is much larger than the 
wire diameter, the phase delay can be approximated as 
proportional to the path length difference (Chuss et al. 
2012 ). 

Ait 7 

0 ( 2 :) — cos6> (3) 
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where A is the wavelength of the radiation and 0 is its 
incident angle on the VPM. 

As can be seen from Equation 2, the VPM modulates 
one of the linear Stokes parameters (in this case, Q) in 
addition to circular polarization (Stokes V). This Q 
V modulation is due to the 45° orientation of the grid 
wires with respect to the detectors. If the grid wires are 
instead at 0° or 90° in this coordinate system, the VPM 
would produce U ^ V modulation, and has no influence 
of the Q signal measured by the signal difference of an 
orthogonal pair of linearly polarized sensors. 

By using sky and/or instrument rotation, the modula¬ 
tion of Stokes Q in the instrument frame can be used to 
measure both Stokes Q and U on the sky. The tracking 
of instrument rotation and pointing position of the in¬ 
strument will also be utilized as secondary modulation 
schemes to track and remove linear polarization from 
ground spill and instrument emission. Because the VPM 
transforms circular polarization to linear polarization, 
VPM modulation can be used to construct a map of 
Stokes V in addition to /, Q, and U. For the CMB, 
the theoretical expectation is that the circular polariza¬ 
tion will be negligible compared to the linear polarization 
(Kamionkowski et al. 1997; Cooray et al. 2003; Alexander 
et al. 2009). There will be additional circular-polarized 
contributions from the atmosphere (Keating et al. 1998; 
Hanany & Rosenkranz 2003; Spinelli et al. 2011; Errard 
et al. 2015) that are not spatially uniform. As such, the 
Stokes V channel will likely be utilized to track these spu¬ 
rious systematic effects, and the sensitivity of the archi¬ 
tecture to astronomical circular polarization will depend 
on the magnitude of such effects. 


2.2. Transmission-Line Model 

To capture details of the VPM operation beyond the 
ideal model, we implement a transmission-line model 
for the VPM based on Chuss et al. (2012). The two 
main manifestations of the improved fidelity of the 
transmission-line model over the ideal model are: (1) 
a small modification of the distance/phase relation in 
Equation 3; and (2) precise treatments of the polariza¬ 
tion isolation and finite absorptance of the VPM’s com¬ 
ponents. This enables the estimation of the differential 
emissivity between the polarization parallel and perpen¬ 
dicular to the wires, which is a key source of systematic 
effects studied here. 


The responses of the wire grid polarizer to parallel and 
perpendicularly polarized incident radiation are modeled 
via inductive and capacitive circuits, respectively. Loss 
in the parallel component is included via a series resis¬ 
tance that is calculable from the wire geometry and mate¬ 
rial. Loss in the perpendicular component can be treated 
in an analogous way, but because it is sub-dominant it is 
not included in this model. 

A combination of the finite conductivity of the wires 
and the imperfect isolation between the reflected and 
transmitted polarizations leads to a higher emissivity for 
the parallel than the perpendicular polarization. In ad¬ 
dition to emission due to its finite physical temperature, 
the VPM also absorbs a small fraction of the signal from 
the sky. The overall instrumental signal is proportional 
to the differential emissivity of the wires and the product 
of the temperature difference between the physical tem¬ 
perature of the VPM and the brightness temperature of 
the background, which in this case consists of contribu¬ 
tions from the atmosphere, the CMB, and foreground 
sources. It can be written as 


sypM = 


TypM 


-( 


^atmTatm T F^nib T ^ ^ 


eVPM 


( 4 ) 

Here, TypM, Fatm, and Tcmb are the physical tempera¬ 
tures of the VPM, atmosphere and CMB, respectively. 
Ti and are the temperature and emissivity of any as- 
trophysical foregrounds in the instrument’s beam. eypM 
is the differential emissivity of the VPM. The result, 
sypM, is the systematic signal emitted from the VPM. 
In this work, we consider the CMB as the only celes¬ 
tial source, so the terms in the summation are set to 
be zero. The magnitude of this systematic signal is a 
function of grid-mirror position. Specifically, the differ¬ 
ential emissivity component, evPM(^) (and consequently 
the systematic signal) is a function of grid-mirror posi- 
tion. For a warm VPM (Tvpm > (eatm^atm + Pcmb)), 
this will lead to a modulated grid emission. For a cold 
VPM (Tvpm < (eatmTatm + Tcmb)), this will lead to a 
modulated grid absorption. 

Selected band-averaged Mueller matrix elements for 
the VPM transmission line model that has been imple¬ 
mented for this work are shown in Figure 3 superposed 
on the corresponding elements for the ideal model. We 
have implemented the transmission line model described 
by Chuss et al. (2012) using a grid having 50 jam di¬ 
ameter wires on a 110 yum center-to-center wire spacing. 
The surface of the wire is copper with a conductivity of 
6 X 10^ S • m“^. Differential emissivity is shown in Fig¬ 
ure 4. A resonant feature in the differential emissivity 
occurs due to the presence of a cavity resonance defined 
by the grid-mirror separation. This resonance is asso¬ 
ciated with the mode having its polarization parallel to 
the wires being trapped by the inductive grid. 

The broadband spectral response of the detectors 
through the VPM is assumed to be a top hat function. 
The frequencies of operation are taken to be between 33 
and 43 GHz (9 to 7 mm), corresponding to the low fre¬ 
quency channel of CLASS (Eimer et al. 2012; Essinger- 
Hileman et al. 2014). In the simulation, the frequencies 
are averaged with a spacing of 0.002 GHz. The feature 
observed in Figure 4 at ^ 3.7 mm is associated with the 
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Fig. 3. — Mueller matrix terms for a VPM described by a 
transmission-line (solid lines) and an ideal (dashed lines) model 
bandpass averaged from 33 to 43 GHz with the wire grid at 45° 
with respect to the polarization sensitive direction of the detec¬ 
tors. The blue and green lines are the Mqq terms, while the red 
and magenta lines are the Mqv terms. The decrease in amplitude 
at larger grid-mirror separations results from the decoherence in 
the modulation across the frequency band of the detectors. 



Fig. 4.— Differential emissivity of the VPM as a function of grid- 
mirror separation band averaged over 33 to 43 GHz. This signal 
leaks into the local Q and U depending on the orientation of the 
grid. The blue line assumes either the detector time constant, Tdet? 
or, equivalently, the modulation frequency of the VPM, /vpm, is 
zero. The red line is the effective (smoothed) differential emissivity 
for the case in which the time constant of the detectors is 0.01 s 
and the modulation frequency of the VPM is 10.1 Hz. This is 
representative of a practical system implementation. 

resonant VPM response in the detector band. As veri¬ 
fication of this physical interpretation, we find this res¬ 
onance disappears when the inductive cavity is removed 
and the polarization parallel to the wires is no longer 
trapped. We verified that the frequency resolution stated 
above was sufficient to Nyquist sample the VPM response 
at all relevant time scales. For a real system, the band¬ 
pass and transfer function of the VPM will need to be 
characterized, but a top-hat frequency response is suf¬ 
ficient to capture the critical aspects of the system re¬ 
sponse. 


Fig. 5.— Difference in the Mueller matrix terms between the 
cases for the grid-detector angle of 45° and 44.5°. Results are for 
Mqq (blue), Mqu (green), Mqv (red). 

3. SYSTEMATIGS 

VPM-based effects, such as VPM emission, are natu¬ 
rally separated from the signal in each detector in the 
following way. As the grid-mirror separation is varied, 
the VPM will differentially affect the linear polarizations 
parallel and perpendicular to its wires. Because the wires 
are at a 45° angle to the polarization sensitive direction 
of the detectors, such linear polarized signals will only 
be present in the local Stokes U parameter, while each 
detector is only sensitive to the local Stokes Q parame¬ 
ter. The CMB signal, on the other hand, will be present 
in the local Stokes /, Q and U parameters. In this sec¬ 
tion, we describe how deviations from this symmetry can 
couple artifacts into the signal. 

3.1. Grid Misalignment Error 

We first examine the effect of a grid misalignment er¬ 
ror. Misalignment arises through an error in the assumed 
angle between the VPM grid wires and the polarization 
axis of the detectors. 

The Mueller matrix of the system is dependent on the 
grid-detector angle. For example, with a perfect 45° 
alignment, Mqq and Mqy are non-zero while Mqu is 
zero. As the angle alignment changes, Mqu becomes 
non-zero, but is still small compared to Mqq and Mqy 
for small deviations from 45°. This is illustrated in Fig¬ 
ure 5 which shows the difference in Mueller matrix terms 
between grid-detector angles of 45° and 44.5°. We test 
the sensitivty of the maps and associated power spectra 
to grid misalignment error. 

3.2. VPM Emission and Absorption 

The main set of systematics that we address in this 
work are those arising from the coupling of the VPM 
emission into the signal. As described in Section 2.2, the 
transmission-line model for the VPM includes the effect 
of the differential emissivity of the wires, which can con¬ 
vert an unpolarized signal into a polarized signal. Since 
the light is emitted, absorbed, and polarized preferen¬ 
tially along the direction of the grid wires, these system¬ 
atics are balanced by placing the grid wires at 45° with 
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respect to the polarization direction of the detectors. In 
this case, in the coordinate system defined by the polar¬ 
ization directions of the detectors (the detector frame), 
the VPM emission and absorption would only affect the 
non-measured U Stokes parameter (i.e. the Mueller ma¬ 
trix term describing / to Q leakage. Mg/, would be zero). 
Under the condition of perfect 45° alignment between the 
grid wires and the VPM, the variation in the I signal in 
each detector pair is common mode. In contrast, the Q 
signal is differential in the two detectors. 

There will always be uncertainty in the orientations 
of the grid wires, and thus they will not be at exactly 
45° with respect to the detectors. If this relative angle 
differs from 45° by an angle 6, Mqj would be non-zero 
and there will be a VPM emission signal of 

Sve = sypM sin S (5) 

leaking to Q in the detector frame, where 5 vpm is the 
total amplitude of the systematic signal emitted by the 
VPM and given by Equation 4. The physical tempera¬ 
tures can be inferred from the data, while eypui^)^ and 
6 are slowly varying and need to be measured. With a 
differential emissivity of 1.0 x 10“^, a temperature dif¬ 
ference of 250 K (Tvpm = 280 K and Tgky = 30 K), and 
6 = 0.05°, the output Q in the detector frame due to 
VPM emission is 2.2 x 10“^ K. This is close to two or¬ 
ders of magnitude larger than the RMS of the Q and U 
CMB polarization signal of ~ 3.6 x 10“^ K, so care must 
be taken in separating the synchronous signature from 
the true sky signal. We discuss this in more detail in 
Section 4. 

As seen in Equation 5, the VPM emissivity amplitude 
is proportional to the difference between the background 
brightness temperature and the physical temperature of 
the VPM. Cooling the VPM to match the brightness tem¬ 
perature of the atmosphere will greatly reduce the ampli¬ 
tude of the VPM emissivity; however, there is a trade-off 
in doing so that is similar to that of other modulators 
(e.g. HWPs). Eor a ground-based implementation, a 
cold modulator is required to be placed behind a pres¬ 
sure window, and thus, polarization effects of the window 
must be considered. Eor space-borne instruments (and 
suborbital platforms such as PIPER), the pressure win¬ 
dow requirement can be relaxed and a cold modulator is 
more pratical. 

In a real VPM-based instrument, the temperature dif¬ 
ference in Equation 5 will vary over time scales that are 
long compared to that of a single stroke of the VPM. 
If the grid is not at 45° with respect to the polarization 
sensitive direction of the detectors, this temperature vari¬ 
ation can introduce noise in the polarized power spectra. 

4. SIMULATIONS 

To investigate the mitigation of the two systematics de¬ 
scribed above, we constructed a simulation of the instru¬ 
ment response during observation. We simulate a CMB 
sky with the parameters from a model consistent with a 
Planck ACDM CMB model, assuming r = 0 (Planck Col¬ 
laboration et al. 2015). We generate a simulated TOD by 
smoothing the simulated sky with a 1.5° Gaussian beam 
and implementing an observing strategy. Next, we gen¬ 
erate maps from these TOD. Einally, we construct the 
power spectra of the maps and subtract the power spec¬ 
tra of simulations without instrumental effects to obtain 


an estimate of the systematic error power spectra. 

The only sources of noise that enter the simulation 
are those from reconstructing the template and from the 
temporal variations in VPM emission signature (see Sec¬ 
tion 5.2 and Section 5.3). Similarly, we ignore the sun, 
foregrounds, and any loss of observation time due to tele¬ 
scope maintenance or calibration to simplify the calcula¬ 
tions. 

To determine the effect of each systematic error term 
on the resulting power spectra, we run several simula¬ 
tions with these effects present. As a representative test 
case, the simulations are modeled on the CLASS experi¬ 
ment and observing strategy. The focal plane simulation 
consists of 36 pixels. Each pixel has two detectors with 
orthogonal polarization orientations, which we designate 
H and V based on their polarization sensitive direction. 
A map of the detector pointing offsets from the boresight 
center of the telescope and the ideal angle of the grid is 
shown in Eigure 6. (This is similar to the 38 GHz focal 
plane for CLASS, though the detectors and VPM polar¬ 
izer in CLASS will be rotated by 45° with respect to the 
directions shown here.) The response of each detector is 
modeled as a single-pole filter with an integration time 
of 10 ms. We scan 360° in azimuth at 40° elevation for 
24 hours. This generates an observed circle whose center 
moves in right ascension over the course of the day. The 
velocity of the scan is varied between 0.5 and 2degs“^ 
the lower and upper velocity limits for CLASS depend¬ 
ing on the declination to make the resulting integration 
time per equal-area pixel as uniform as possible over the 
CLASS survey area. Each simulation consists of 7 days 
of observations. The boresight angle of the telescope is 
rotated by 15° every day, through a total of 90° of rota¬ 
tion. This will reduce systematics due to polarized emis¬ 
sion that is fixed in the instrument frame. To be conser¬ 
vative, we chose these values so that, due to symmetry, 
the systematic error is maximized in the reconstructed 
U maps and canceled out in the reconstructed Q maps. 
This maximizes the effect of the systematic in the BB 
power spectrum. Our coordinate system is chosen such 
that at a boresight angle of 0° and a parallactic angle of 
0°, Q in the detector frame is the same as Q in the sky 
frame. 

Eor this paper, we first difference the H and V TOD 
before making maps. After differencing the TOD, we 
then solve the standard least-squares map-making equa¬ 
tion (Tegmark 1997): 

(A^N-iA)x = A^N-M (6) 

where N is the noise covariance matrix, d are the mea¬ 
sured TOD, A is the pointing and modulation encoding 
matrix, and x is a vector containing the Stokes parame¬ 
ters in each map pixel. 

The rapid polarization modulation significantly re¬ 
duces the cross-coupling between pixels on the sky com¬ 
pared to scanning techniques. Euture work will quanti¬ 
tatively address residual correlations; however, to target 
the specific systematic effects of interest here, we consider 
the noise to be constant and uncorrelated, allowing a use¬ 
ful simplication of Equation 6. If N is a constant multi¬ 
plied by the identity matrix then (A^N“^A) is a block 
diagonal matrix where there are no correlations between 
map pixels. The Stokes parameters in each map pixel can 
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Fig. 6.— The focal plane used in the simulations. The lines show 
the polarization directions of each detector used in the simulations. 
There are 36 pixels containing 72 polarization sensitive detectors. 
The lines in the upper left corner show the ideal direction that the 
grid makes relative to the polarization sensitive direction of the 
detectors. 


Fig. 7. — The amplitude of the VPM modulation of the polariza¬ 
tion signal as a function of time. The black dots are the positions 
sampled in the simulation due to the 100 Hz sampling rate over the 
first 0.3 s. The VPM is modulated at a frequency of 10.1 Hz. The 
VPM stroke profile is chosen to maximize sensitivity to Stokes Q 
at the expense of sensitivity to Stokes V. 


therefore be solved for independently of all other pixels. 

The pointing matrix, A, encodes the pointing along 
with how the associated signal transforms as a function 
of the parallactic angle and VPM position. For an ex¬ 
periment with a VPM, the simplest model of the TOD, 
when the grid-detector angle is at 45°, is 

d{z, a) = QMqq cos 2a+ U Mqq sin 2a+ V Mqv (7) 

where the QQ, and QV elements of the VPM Mueller 
matrix (M) are dependent on the grid-mirror separa¬ 
tion, 2 ;. The parameter a is the sum of the parallactic 
angle and Foresight angle of the telescope and is the an¬ 
gle between the detector coordinate system and the sky 
coordinate system. A more accurate model for the signal 
that is valid for general grid-detector angles is 

d = QMqq cos 2a - QMqu sin 2a+ U Mqq sin 2a 
+ UMqu cos 2a + V Mqv (8) 

and incorporates the Mqu VPM Mueller matrix term 
that is non-zero when the grid-detector angle is not at 
45°. Plots of Mqq and Mqv for a grid-detector angle 
of 45° as a function of grid-mirror separation are shown 
in Figure 3. The grid-mirror separation as a function of 
time is given by z{t) = 2^0 + ^ [1 — cos(27r/t)], where 
zq = 1850 jam is the minimum grid-mirror separation, 
Az = 2640/im is the mirror throw, and / = 10.1 Hz 
is the frequency of the mirror modulation rate. The 
grid-mirror separation range is chosen to maximize sig¬ 
nal to noise, and the modulation frequency is chosen to 
be slightly different from 10 Hz to provide more than 10 
unique positions given the 100 Hz sampling rate of the 
data. With this throw, we include the resonant feature 
in the smoothed (red) curve for differential emissivity 
visible in Figure 4. The corresponding transfer functions 
as a function of time, MQQ(t) and MQv(t), are shown in 
Figure 7. The functions shown here are calculated using 
the transmission-line model of the VPM. In practice, the 
fidelity of a real experiment can be improved through 
measurement of the VPM transfer function. The point¬ 


ing matrix. A, encapsulates the coefficients of the three 
polarization signals: Q, V, and V. This technique was 
used in the making of Figure 1. 

After making maps, we construct the auto-correlation 
power spectra using standard HEALPix^ routines. 
While we measure nearly the full sky, there is some leak¬ 
age from E to B due to the sky cut if the standard 
pseudo-C^ estimator is used. We implement the pure 
Ci estimator (Smith 2006; Grain et al. 2009) to mini¬ 
mize this leakage. The residual leakage is several orders 
of magnitude below r = 0.01 in the simulations so any 
significant level of B due to systematics is not associ¬ 
ated with this reconstruction. The circular polarization 

V is treated the same way as the temperature anisotropy 
when constructing the power spectra, because while a 
rotation of the coordinate system will mix Q and V, it 
does not affect the measured value for either the / or 

V Stokes parameters. Therefore the TT and VV power 
spectra will have similar mode mixing properties when 
examining maps having partial sky coverage. 

The C£ power spectra shown throughout this paper are 
the difference between the reconstructed power spectra 
with and without systematics averaged over 100 realiza¬ 
tions. This comparison is done in power spectrum space 
to provide sign dependence for the E-mode systemat¬ 
ics. By differencing power spectra, we can distinguish 
between systematics that increase the overall power and 
systematics that might shift power from E-modes to B- 
modes and therefore decrease the overall amplitude of 
the EE power spectrum. 

While we show results from simulations with no input 
E-modes (r = 0), all reconstruction techniques are also 
applied to simulations with a non-zero input BB power 
spectrum and no systematics. We have found that the 
reconstruction techniques studied in this paper do not 
bias the underlying measurements, as we find that the 
resulting BB power spectrum is unbiased compared to 
the input BB power spectrum. The only bias is due to 

^ http://healpix.sf.net 








the partial sky cut and that is removed by the imple¬ 
mentation of the pure Cg estimator and comparisons to 
simulations with no systematics. The techniques do not 
remove any measurable fraction of the input 5-modes or 
convert input 5-modes to 5-modes. 

5. RESULTS 

We add in the systematics in three cumulative steps. 
First, we include the effect of grid misalignment without 
the effect of VPM emission. After understanding these 
results, we add in a non-zero VPM emission. Finally, in 
the last simulation we allow the temperature difference 
between the VPM and the atmosphere to vary over time. 

5.1. Grid Misalignment Error 

Reconstructed Q, 5, and V maps for grid misalign¬ 
ment error are shown in Figure 8. The grid is set to 44.5° 
when simulating the data and assumed to be 45° when 
reconstructing the maps. Here, we do not include the 
effects of differential emissivity (ie., the VPM physical 
temperature is set to the brightness temperature of the 
sky in Equation 5). An error of 0.5° is a conservative es¬ 
timate for the uncertainty in grid-detector alignment as 
this angle will be measurable to greater precision than 
this. There is no visible systematic error dominating the 
reconstructed Q and U maps; however, the V map shows 
leakage from the input Q map, with the sign depending 
on the declination. This is due to the slight difference 
in the V modulation term with the grid at 44.5° and 
45.0°. The sign is linked to the scan strategy and how 
the orientation of the polarization axis of the detectors 
with respect to the sky coordinate system changes over 
the course of the day. 

The amount of Q and U mixing is too small to be 
evident in just the reconstructed maps; however, some 
mixing does occur and is evident in the difference maps. 
The AQ map looks like the U map and the A5 map looks 
like —Q. Because of the misalignment, the measurement 
of the local Q parameter includes a small contribution 
from the local U due to the extra VPM modulation of 
U into Q. Due to the parallactic angle variation this 
results in Q in the sky coordinate system having some 
contribution from U and vice-versa. In addition, there 
is a residual coupling into the reconstructed V map due 
to small differences between the expected and actual V 
modulation. 

The resulting systematic power spectra are shown in 
Figure 9. The systematic 55 power spectrum is ~ 10% 
of the 5-mode power for r = 0.01 on the largest scales. 
The average systematic over the range 10 < ^ < 100 is 
a few percent of r = 0.01, though it does peak above 
r = 0.01 for 100 < ^ < 200. Both the 55 and VV 
systematic power spectra have shapes similar to the EE 
primordial power spectrum and we can be certain that 
the 5-modes are leaking to 5 and V due to the grid 
misalignment. In addition a small amount of 55 would 
be leaking to EE and W, but its fractional amplitude 
is small and is not an issue. This is expected since the 
mixing of Q and U in the maps is similar to the case of 
polarization rotation. The systematic EE spectrum is 
negative because we are comparing the recovered power 
spectrum to one calculated in a simulation without sys¬ 
tematics. Since some E leaks to 5 and V, the EE power 


spectrum exhibits a decrease in amplitude relative to the 
case in which the grid misalignment is absent. 

The 5 to 5 leakage introduced by the grid mis¬ 
alignment will generate non-zero TB and 55 cross¬ 
correlation power spectra. This systematic effect can be 
removed by either implementing a self-calibration tech¬ 
nique that employs an in-situ measurement of the ro¬ 
tation angle using the cross-correlation power spectra 
(Keating et al. 2013) or by directly correcting the transfer 
function given a precise measurement of the grid-detector 
angle. This systematic is due to the fixed misalignment 
between the actual grid-detector angle and the value we 
used in the transfer functions in the simulation. This 
single degree of freedom can be measured and corrected. 
For small errors in the grid-detector alignment, the leak¬ 
age of 5 to 5 depends linearly on the size of the mis¬ 
alignment. It should be straightforward to correct for the 
error and lower the amplitude of this systematic below 
the amplitude of the noise. 

5.2. VPM Emission 

Our approach for mitigating VPM emission is to re¬ 
move its signature from the TOD. To do this, it is neces¬ 
sary to understand the functional form of the emission as 
a function of grid-mirror separation. In this section, we 
describe the generation of the VPM emission template 
and then show that it can be used to clean the TOD in 
the case of a constant environment. In the next section, 
we consider VPM emission removal in the presence of 
varying VPM/atmospheric temperature. 

One method of template reconstruction is to lock the 
detector sampling rate to the VPM throw (e.g. /vpm = 
10 Hz, /d = 100 Hz). In this case, a set number of discrete 
value of the grid-mirror separation would be repeatedly 
sampled. Because we are using the actual data to extract 
the VPM emission template, these would be the same 
positions at which the removal of the template would 
ultimately take place. 

To simulate the effect of VPM emission, we set the an¬ 
gle between the VPM and the detectors to be 44.5°, as in 
the previous section, and set the temperature difference 
between the VPM and the sky to a constant 280 K. 

Reconstructed maps are shown in the top row of Fig¬ 
ure 10. To remove the VPM emission from the recon¬ 
structed maps, we construct a template of the VPM emis¬ 
sion using the TOD, and use this template to remove the 
VPM emission from the TOD. We construct the tem¬ 
plate by partitioning the grid-mirror separation into a 
number of bins and determining the average signal over 
all detectors for the entire observation in each bin. The 
sky signal is dependent on the position of the detectors 
and the polarization orientation on the sky, whereas the 
VPM emission is only a function of grid-mirror separa¬ 
tion. Therefore, averaging measurements over the entire 
observation results in the sky signal averaging to zero 
with the residuals being a good estimate of the VPM 
emission template. This template is then interpolated to 
the observed grid-mirror positions using a cubic interpo¬ 
lation. For grid-mirror positions that are larger than the 
average value of the last bin or smaller than the average 
value of the first bin and therefore outside the range of in¬ 
terpolation, we apply a quadratic extrapolation. Finally, 
we subtract the interpolated template from the TOD of 
each detector. 
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Fig. 8.— The top row is the input Q, t/, and V maps for one realization. The middle row is the reconstructed Q, U^ and V maps for a 
single realization of the simulation with the VPM grid at 44.5° and no VPM emission. The bottom row is the difference between the top 
two rows. 



Fig. 9. — Systematic power spectra for the grid misalignment 
error for a grid-detector angle of 44.5°. The red line is the EE 
systematic power spectrum. The green line is the BB systematic 
power spectrum. The blue line is the VV systematic power spec¬ 
trum. The solid lines are positive and the dashed lines are negative. 
The black lines are EE (upper) and BB (lower) power spectra for a 
ACDM model with r = 0.01. The systematic power spectra are the 
difference between the power spectra constructed from simulations 
with the systematic present and ones without any systematics. 


Reconstructed maps with the template subtraction are 
shown in the bottom row of Figure 10. The template used 
in this simulation has 2000 bins, 500 of which contain 
more than one measurement in the 24 hour simulation. 
This number is realizable given the > 2600 /im throw and 
a possible 0.1 /im resolution of the VPM grid-mirror sep¬ 
aration. We find that the number of bins can be reduced 
by close to a factor of 10 while still maintaining the abil¬ 
ity to remove the systematic. Plots of the reconstructed 
power spectra in Figure 11 demonstrate that the tem¬ 
plate can be used to successfully reconstruct and remove 
the VPM emission. The dominant remaining systematic 
in EE and BB power spectra is the leakage due to the 
grid misalignment. A comparison with Figure 9 illus¬ 
trates that any residual VPM emission is significantly 
below the r = 0.01 level. 

Though system (detector and photon) noise for the in¬ 
dividual sensors in the instrument is not included in this 
simulation, it is important to check that these factors will 
not significantly limit this technique. First, the noise in¬ 
troduced by the interpolation can be estimated by the 
variance between the modeled VPM emission and the 
reconstructed template, which is found to be ~ 300 nK. 
This is significantly below the anticipated sensor noise, 
200 /iKs^/^, in a typical channel. Thus, we do not expect 
the template reconstruction to significantly decrease the 
sensitivity of the experiment. Second, this technique as¬ 
sumes that the signal-to-noise is large over time scales 
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short compared to changes in the brightness tempera¬ 
ture of the atmosphere. We anticipate a detection of the 
^ 100 /jK VPM emission signal on ^minute time scales. 
This is significantly below the total integration time used 
to calculate the VPM emission, and thus we anticipate 
that the presence of noise will not significantly adversely 
affect the result of the template reconstruction described 
here. This estimate implicitly assumes that the reference 
phase can be derived from the emission signature without 
incurring additional noise on the time scales of interest. 
Alternatively, the modulation reference and data acqui¬ 
sition can be synchronized to eliminated this potential 
concern - in practice this would be the preferred experi¬ 
mental approach. 

5.3. Variable VPM Emission 

The amplitude of the VPM emission is directly pro¬ 
portional to the difference between the temperature of 
the VPM and the brightness temperature of the back¬ 
ground, which is dominated by the atmosphere. There¬ 
fore, the temperature excursions in either component will 
cause the emission signal to vary in time. We generically 
model these time-dependent temperature excursions us¬ 
ing a 1//^ power spectrum. The amplitude is chosen 
to be 2.5 X 10“^K/\/Hz at IHz. This leads to a peak- 
to-peak variation in the VPM temperature difference of 
around 2.5 K when a 1.0 x 10“^ Hz low frequency cutoff 
set by calibration considerations is applied (i.e. some¬ 
thing is needed to regularize the spectral powers). In 
practice, these variations will be dominated by the at¬ 
mospheric contribution, and thus this noise spectrum is 
effectively modeling fluctuations in the brightness tem¬ 
perature of the atmosphere. The VPM fluctuations are 
expected to be much smaller. Components deployed in 
similar thermal environments have observed temperature 
drifts of 0.5 K on time scales of hours under passive ther¬ 
mal control (Wollack et al. 1997). A real implementation 
of a VPM should include temperature sensor data that 
can be used in the analysis to understand and remove 
fluctuations in the VPM emission due to a varying tem¬ 
perature. Such secondary temperature tracking is not 
considered in this analysis, as the focus here is on the 
removal of this systematic using radiometric information 
from the TOD alone. 

With a time-varying temperature difference between 
the VPM and the sky, the constructed template is the 
VPM emission at the mean temperature. At any point in 
time, there will be some residual VPM emission whenever 
the temperature difference is not at its mean value. This 
residual emission may result in an unacceptable level of 
a systematic BB power spectrum and it is necesary to 
correct for this temperature dependence. 

For this simulation, we correct for the grid misalign¬ 
ment in the polarization transfer function by using Equa¬ 
tion 8 as the model for the TOD. A plot of the recon¬ 
structed power spectra with mean teperature template 
subtraction is shown in Figure 12. The resulting power 
spectrum is above r = 0.01 between ^ = 10 and = 20 
and would be a concern for measuring r on large angular 
scales. The amplitude will decrease with additional data 
since the temperature is modeled as a 1//^ noise power 
spectrum. 

To correct for the residual 1//^ noise from the temper- 


Sample Rate 

Value 

Data Sample Rate 

100 Hz 

VPM Sample Rate 

10 Hz 

Beam Scan Rate 

0.5 - l.OHz 

Destripe Rate 

0.05 Hz 

Calibration Freqeuency 

10-4 Hz 


TABLE 1 

The sampling rates for the various levels of the 

MODULATION PROCESS. 


ature dependence of the VPM emission, we adopt stan¬ 
dard destriping techniques (Delabrouille 1998; Revenu 
et al. 2000; Kurki-Suonio et al. 2009; Keihanen et al. 
2010) to determine the time-dependent temperature dif¬ 
ference between the VPM and the sky from the TOD. 
In this case, we have two sets of destriping amplitudes 
that must be determined for a single length of time. The 
VPM emission as a function of the VPM temperature 
can be written as a linear combination of a set of ampli¬ 
tudes multiplying the VPM emission at 280 K (template 
amplitudes) and another set of amplitudes that are off¬ 
sets to the TOD (offset amplitudes). In other words, the 
VPM emission at any temperature can be written as a 
function of the VPM emission at 280 K as 

5ve(T) = a(T)5ve(280K) + b{T) (9) 

where a are the template amplitudes and b are the off¬ 
sets. Both are functions of the temperature. The only 
difference between the template amplitudes and offsets 
in the absence of other 1//^ noise will be a multiplica¬ 
tive constant as both amplitudes should reconstruct the 
temperature change over time. Our implementation of 
destriping incorporates these two sets of amplitudes. For 
the purposes of this simulation, we fit for amplitudes in 
20 s intervals. There are « 200 full cycles of the VPM 
over this interval. The sampling rates for the various 
levels of the modulation process are shown in Table 1. 

The full derivation of the destriping procedure is dis¬ 
cussed in depth elsewhere (Kurki-Suonio et al. 2009). We 
summarize the equation for the optimal amplitudes here. 
The two sets of amplitudes are the results of solving the 
equation 

(F^N-iZF)a = F^N-^Zd (10) 

where 

Z = I- 

M = A^N“iA. (11) 

A, d, and N were defined for Equation 6 . I is the 
identity matrix and F is the matrix that expands the 
template and offset amplitudes to the TOD. The vari¬ 
able Z projects out a best estimate of the sky signal. 
The optimal amplitudes, a, are determined by solving 
Equation 10 using a preconditioned conjugate gradient 
method. The best fit VPM emission model is then re¬ 
moved from the TOD before estimating the map as in 
Equation 6. 

Power spectra are shown in Figure 13 for a simula¬ 
tion in which the TOD are corrected using the template 
subtraction combined with the implementation of the de¬ 
striping algorithm described above. In addition, the po¬ 
larization transfer function is corrected for grid misalign¬ 
ment by including the correct Mqq and Mqu Mueller 
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Fig. 10.— Reconstructed Q, U, and V maps for a single realization of the simulation with the VPM grid at 44.5° and VPM emission 
from a VPM at 280 K. The first row has simulations without VPM emission template subtraction from the TOD. Such simulations are 
dominated by the VPM emission systematic. The second row subtracts the VPM emission residuals from the TOD. 




Fig. 11.— Systematic power spectra for VPM emission with 0.5° misalignment and a constant VPM temperature. (Left) We do not 
subtract a template from the TOD. (Right) We subtract the VPM emission template from the TOD. In this case, BB systematics are 
dominated by grid misalignment error as in Figure 9. The systematic power spectra are the difference between the power spectra constructed 
from simulations with the systematic present and ones without any systematics. 


matrix terms in reconstruction. In this case, errors in 
the interpolation of the binned template that couple to 
the time varying temperature necessitate the use of the 
2000 bin template. The resulting BB power spectrum is 
more than an order of magnitude below the target level 
of r = 0.01. On the largest scales {i < 30), the residual 
systematic is on the level of a grid misalignment of 0.1°. 

We also consider the case for which the VPM frequency 
is synchronized to the detector data rate. In this case, 
we take a VPM frequency of 10 Hz, and detector rate 
of 100 Hz, so that the detector is sampling the same 10 
positions of the VPM throw for every stroke. In this 
case, no interpolation is required and the map residuals 
are similar to those achieved with the 2000 bins and an 


arbitrary sampling scheme. 

6. CONCLUSIONS 

We have run simulations to study how systematics and 
noise associated with the thermal variation of the atmo¬ 
sphere and VPM will affect measurements of the CMB 
and have highlighted which aspects of a VPM-based ex¬ 
periment need to be measured and/or controlled to reach 
a desired value of r. We have examined the effect of an er¬ 
ror in the VPM transfer function (through an error in the 
detector-grid alignment), VPM emission and absorption, 
and temperature variation of the VPM and atmosphere 
on the reconstructed Q and U maps and have propa¬ 
gated these results to the BB power spectrum. None of 
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Fig. 12.— Systematic power spectra for a simulation with VPM 
emission and a varying VPM temperature. The systematic power 
spectra are the difference between the power spectra constructed 
from simulations with the systematic present and ones without 
any systematics. An average template as a function of VPM grid- 
mirror separation is subtracted from the TOD to remove the VPM 
emission. Variation in the VPM emission about the average results 
in systematics comparable to 5-modes with r = 0.01. It is possible 
to mitigate this systematic using the destriping algorithm described 
in the text. (See Figure 13.) 



Fig. 13.— Systematic power spectra for a simulation with VPM 
emission and a varying VPM temperature. The systematic power 
spectra are the difference between the power spectra constructed 
from simulations with the systematic present and ones without 
any systematics. The destriping algorithm is used to correct the 
template for the temperature variation to remove the systematic 
from the TOD. The contamination to BB is seen to be an order of 
magnitude below the r = 0.01 level. 


the systematics studied fundamentally limit the detection 
of r on large scales for r = 0.01. 

The simulations show that a systematic due to grid- 
angle error may need to be treated for high-fidelity char¬ 
acterization of inflationary 5-modes. An error of 0.5° 
in the simulation can be 10% of the 5-mode power at 
r = 0.01. Smaller and more realistic angular errors will 
produce insignificant contamination. In the simulations 
with a time varying VPM temperature, we removed the 
grid-angle error systematic by incorporating a non-zero 


Mqu term into the TOD model, removing most of the 
grid-angle error systematic. Measurements of the VPM 
transfer function are needed to ensure that this 5 to 5 
leakage is minimized. Alternatively, the grid-angle sys¬ 
tematic can be determined by fitting for a rotation angle 
from the TB and EB power spectra and correcting the 
transfer function similar to the self-calibration procedure 
of Keating et al. (2013). This would preclude a mea¬ 
surement of isotropic cosmic birefringence, so a hardware 
calibration system (Johnson et al. 2015) would likely be 
implemented. 

We have shown that the temperature to polarization 
leakage due to VPM emission can be mitigated by re¬ 
moving a template from the TOD. The VPM emission 
is a function of grid-mirror separation and a high signal- 
to-noise template can be constructed since we repeatedly 
observe at every different VPM position. Different tem¬ 
plates may be needed for each separate pixel depending 
on how the VPM emission changes in detail across the 
focal plane. The template-based emission model remains 
effective when allowed to vary in time. Implementing de¬ 
striping with a set of amplitudes that multiply the tem¬ 
plates allows us to separate most of the time varying 
systematic from the CMB signal. Additionally, cross¬ 
spectra between different seasons may alleviate the need 
for the destriping as the temperature variation would be 
uncorrelated between the two maps. As long as the av¬ 
erage template is removed the temperature variations 
would only add noise to the power spectra instead of 
extra signal. 

While the results shown in this paper assume that the 
data from the pair of orthogonally-polarized detectors in 
a pixel are differenced, the rapid modulation by the VPM 
will allow polarization maps to be made from the undif¬ 
ferenced data from a single detector and will allow pixels 
with only one working detector to be used in the analy¬ 
sis. The ability to separate systematics from the signal 
in this case is currently being investigated and prelimi¬ 
nary results show that it should not limit the sensitivity 
of the experiment. 

The use of a VPM will enable limits on the Stokes V 
parameter. This is expected to be a null map in standard 
ACDM cosmology, but there may be some foregrounds 
(Keating et al. 1998; Spinelli et al. 2011) that can be 
measured. The introduction of the V Stokes parame¬ 
ter results in four additional power spectra that we can 
construct: VV, TV, EV, BV. We have determined the 
systematic leakage to the VV power spectrum in the sim¬ 
ulations. The amplitude of the systematic is on a level 
similar to that for EE and BB. For the system and envi¬ 
ronmental limits considered here, the limit on detecting 
circularly polarized foregrounds is ^ 100 nK. This is sev¬ 
eral of orders of magnitude better than the current 95% 
upper limits on the degree of the CMB circular polariza¬ 
tion (Mainini et al. 2013). 
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